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Experiments and simulations have been performed to examine the finite-temperature behavior of 
Gai7^ and Ga2o^ clusters. Specific heats and average collision cross sections have been measured as 
a function of temperature, and the results compared to simulations performed using first principles 
Density-Functional Molecular-Dynamics. The experimental results show that while Gair^ appar- 
ently undergoes a solid-liquid transition without a significant peak in the specific-heat, Ga2o''' melts 
with a relatively sharp peak. Our analysis of the computational results indicate a strong correla- 
tion between the ground-state geometry and the finite-temperature behavior of the cluster. If the 
ground-state geometry is symmetric and "ordered" the cluster is found to have a distinct peak in 
the specific-heat. However, if the ground-state geometry is amorphous or "disordered" the cluster 
melts without a peak in the specific-heat. 

PACS numbers: 61.46.-l-w, 36.40.-c, 36.40.Cg, 36.40.Ei 
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I. INTRODUCTION 

It is now well established that the mehing points of 
particles with thousands of atoms decrease smoothly with 
decreasing particle size, due to the increase in the sur- 
face to volume ratioi^ However, unlike particles with 
thousands of atoms or the bulk material, probing the 
finite-temperature properties of small clusters (with < 
500 atoms) is non-trivial and remains a challenging task. 
Experimental studies of the melting transitions of clus- 
ters in small size regime have only recently become pos- 
siblei^*^i^*S*iSi2iiSiii Several interesting phenomena have 
been observed, including melting temperatures that rise 
above the bulk valueS^ and strong size-dependent vari- 
ations in the melting temperatures ."^^9. These exper- 
imental findings have motivated many theoretical in- 
vestigations on the finite-temperature behavior of clus- 
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Simulations based on first principles 
have been particularly successful in quantitatively ex- 
plaining the factors behind the size-dependent variations 
in the melting behavior of clusters lii Thus, a confluence 
of recent advances in experimental methods and theoret- 
ical studies using first principles methods have set the 
stage for a major increase in our understanding of phase 
transitions in these small systems. 

Coming to the present work on gallium clusters, it 
is by now well known that gallium clusters not only 
melt at substantially higher temperatures than the bulk 
{Tjnihuik] = 303 K),- but they also exhibit wide variations 
in the temperature dependences of their specific-heats, 
with some clusters showing strong peaks (due to the la- 
tent heat), while others (apparent " non-melters" ) show- 
ing no peaki^i. These features show a strong dependence 
on cluster size, where the addition of a single atom can 
change a cluster with no peak in the specific-heat into a 



"magic melter" with a very distinct peak. This behavior 
has been observed for gallium clusters, Ga„, with n — 
30-55. 

In the present work, we probe the melting behavior 
of small gallium cluster ions where we show that the 
"non-melting" and "melting" features in the specific- 
heats are observed in clusters as small as Gan~^ and 
Ga2o~'', respectively. Prior experimental results for Ga]^y 
over a limited temperature range showed no evidence for 
a melting transition^ The experimental results in this 
case were specific-heat measurements performed using 
multi-collision induced dissociation, where a peak in the 
specific-heat due to the latent heat was the signature of 
melting. On the other hand, recent simulations for Gai7 
show a broad peak in the specific-heat centered around 
600 K. The previous specific-heat measurements for Ga^y 
extended only up to 700 K, so one possible explanation 
for this apparent discrepancy is that the melting transi- 
tion occurred at a slightly higher temperature than ex- 
amined in the experiments. Here, we report specific-heat 
measurements for Ga^j^y over a more extended tempera- 
ture range, along with specific-heat measurements for 
Ga^Q. While no peak is observed in the heat-capacities 
for Ga^y, a peak is observed for Ga^Q. 

To further probe the melting transitions in these clus- 
ters, ion mobility measurements were performed for Gaj^y 
and Gajg as a function of temperature. The ion mobil- 
ity measurements provide average collision cross sections 
which can reveal information about the shape and vol- 
ume changes that occur on melting. For example, a clus- 
ter with a non-spherical geometry might be expected to 
adopt a spherical shape (a liquid droplet) on melting. If 
there is not a significant shape change, there may still 
be a volume change on melting. Most bulk materials ex- 
pand when they melt (the liquid is less dense than the 
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solid). Even in the absence of a significant shape or vol- 
ume change, the cross sections might show an inflection 
at the melting transition due to the thermal coefhcient of 
expansion of the liquid cluster being larger than for the 
solid (in the macroscopic regime most liquids have larger 
coefficients of expansion than the corresponding solids). 
An inflection is observed in the cross sections for both 
Gsl^-j and Ga^Q. Thus, the ion mobility measurements 
suggest that Ga^t, as well as Ga^Q are in a liquidlike state 
above 800 K. 

To explore the reasons behind the behavior out- 
lined above (i.e., to determine why Ga^^ apparently 
melts without a peak in its specific-heat, while a 
peak is observed for Ga^p), we have carried out first 
principles Density-Functional (DF) Molecular-Dynamics 
(MD) calculations on both clusters. The ground-state 
structure and the bonding within the clusters is ana- 
lyzed. The ionic specific-heat is computed using multiple 
histogram method. The calculated specific-heats for 
Ga^7 show three broad low intensity maxima that ex- 
tend from 300 to 1400 K. This resembles the experiment 
results where the measured specific-heats are relatively 
featureless. In contrast, the calculated specific-heats for 
GaJo show a clear peak around 750 K. This is in excel- 
lent agreement with the peak obtained from experimen- 
tal measurements (which occurs at around 700 K). Fi- 
nally, our theoretical results show that the features in the 
specific-heat curves are influenced by the ground-state 
geometry, the bonding of the atoms within the ground- 
state structure, and the isomer distribution that becomes 
accessible as the temperature is raised. 

In Sec. ^ we present the experimental methods and 
the computational details. In Sec. Illll we discuss the ex- 
perimental and theoretical results on both clusters. We 
conclude our results in Sec. IIVI 



II. METHODOLOGY 

Specific-heats were measured using the recently devel- 
oped multi-collision induced dissociation approach. The 
cluster ions are generated by laser vaporization of a liquid 
gallium target in a continuous flow of helium buffer gas. 
After exiting the laser vaporization region of the source 
the clusters travel through a 10 cm long temperature vari- 
able extension where their temperature is set. Cluster 
ions that exit the extension are focused into a quadrupole 
mass spectrometer where a particular cluster size is se- 
lected. The size selected clusters are then focused into a 
collision cell containing 1 Torr of helium. As the clusters 
enter the collision cell they undergo numerous collisions 
with the helium, each one converting a small fraction of 
the ions translational energy into internal energy. If the 
initial kinetic energy is high enough some of the cluster 
ions may be heated to the point where they dissociate. 
The dissociated and undissociated cluster ions are swept 
across the collision cell by a small electric fleld and some 
of them exit through a small aperture. The ions that exit 



are analyzed in a second quadrupole mass spectrometer 
and then detected by an off-axis collision dynode and 
dual microchannel plates. The fraction of the ions that 
dissociate is determined from the mass spectrum. Mea- 
surements are performed as a function of the ions initial 
kinetic energy, and the initial kinetic energy required for 
50% dissociation (IKE50%D) is determined from a lin- 
ear regression. IKE50%D is measured as a function of 
the temperature of the temperature-variable extension 
on the source. IKE50%D decreases as the temperature 
is raised because hotter clusters have more internal en- 
ergy, and hence less energy needs to be added in order 
to cause dissociation. At the melting transition a sharp 
decrease in IKE50%D is expected due to the latent heat. 
The derivative of IKE50%D with respect to temperature 
is approximately proportional to the specific-heat. The 
proportionality constant is the fraction of the clusters ini- 
tial kinetic energy that is converted into internal energy, 
which is estimated from an impulsive collision model. A 
drop in the IKE50%D values due to the latent heat of a 
melting transition leads to a peak in the specific-heat. 



Ion mobility measurements can provide information on 
the shape and volume changes that occur when clusters 
melt. For the ion mobility measurements, the collision 
cell is replaced by a 7.6 cm long drift tube. 50 fis pulses 
of cluster ions are injected into the drift tube and the drift 
time distribution is obtained by recording the ions arrival 
times at the detector with a multichannel scalar. Average 
collision cross sections are obtained from the drift time 
distributions using standard mcthodsi^ 



All the simulations are performed using Born- 
Oppenheimer molecular-dynamics based on Kohn-Sham 
formulation of Density-Functional Theory (DFT)i2S 
We have used Vanderbilts' ultrasoft pseudopotentials^ 
within the GGA approximation, as implemented in the 
VASP package24 for both clusters. For all calculations, 
we use only 4s^ and 4p^-electrons as valence, taking the 
3ci-electrons2^ as a part of the ionic core. An energy 
cutoff of about w 10 Ry is used for the plane-wave ex- 
pansion of the wavefunction, with a convergence in the 
total energy of the order of 0.0001 eV. Cubic supercells 
of lengths 20 and 25 A are used for Gaf^ and Ga^g, re- 
spectively. For examining the finite-temperature behav- 
ior, the ionic phase space of the clusters is sampled by 
isokinetic MD where kinetic energy is held constant via 
a velocity scaling method. For both the clusters, we split 
the total temperature range from 100-1400 K into 15 
different temperatures. We maintain the cluster at each 
temperature for a period of at least 90 ps, leading to 
total simulation times of the order of 1 ns. The result- 
ing trajectory data were used to compute standard ther- 
modynamic indicators as well as the ionic specific-heat, 
via a multihistogram technique. Details can be found in 
Bel. 120112ft 



III. RESULTS AND DISCUSSION 



Specific-heats measured for Ga^y and Ga^Q as a func- 
tion of temperature are sliown in the lower half of Fig. ^ 
The points are the experimental values, while the dashec 
line is the prediction of a modified Debye model. In the 
case of Ga^7, the specific-heats shown in Fig. Qappeai 
to gradually increase up to around 900 K. The sharp de 
crease in the specific-heats above 900 K is an artifaci 
due to evaporative cooling, the spontaneous unimolec 
ular dissociation of the cluster ions as they travel be- 
tween the source extension and the collision cell. Foi 
Ga^Q, the specific-heats show a broad maximum, arouna 
400 K wide, centered at around 725 K. The peak for Gajg 
is significantly broader than observed for larger clusters 
(like Gajj and Ga^g) where the peak was attributed to 
a melting transition. However, it is well known that the 
melting transition, and the corresponding peak in the 
specific-heats, becomes broader with decreasing cluster 
size. Thus, even though the peak in the specific-heats 
for Ga^Q is around 400 K wide, it is appropriate to as- 
sign it to a finite-size analog of a bulk melting transition. 
The center of the peak is at around 725 K, this is well 
above the bulk melting point (303 K). This continues a 
trend reported for larger cluster sizes (n = 30-55) where 
the melting temperatures are also significantly above the 
bulk value. The unfilled red circles in Fig. show the 
average collision cross sections determined for Gaf^ and 
Ga^Q as a function of the temperature. The cross sections 
are expected to systematically decrease with increasing 
temperature because the long range attractive interac- 
tions between the cluster ion and the buffer gas atoms 
becomes less important, and the collisions become harder 
as the temperature is raised. The thick dashed red line 
in the figures show the expected exponential decrease in 
the cross sections with increasing temperature. There is 
an infiection in the cross sections for Gajg that appears 
to slightly precede the peak in the specific-heat for this 
cluster. The inflection in consistent with a melting tran- 
sition where the liquid cluster has a larger coefficient of 
thermal expansion than the solid. There is also an infiec- 
tion in the cross sections for Ga]'^. This suggests that a 
solid-liquid transition also occurs for Gafy, but without 
a significant peak in the specific-heat. 

To understand the reason behind the different behavior 
observed for Gair"*" and Ga2o''', we have carried out a de- 
tailed analysis of structure and bonding in both clusters. 
As will be become apparent from the following discus- 
sion, that the ground-state geometry and the nature of 
bonding plays a crucial role in determining the finite- 
temperature behavior of the cluster. We begin with 
a discussion of the ground-state geometries of cationic 
Gai7 and Ga2o clusters. We have obtained more than 50 
distinct equilibrium geometries by quenching more than 
200 structures, selected from a few high temperature MD 
runs, for both sizes. In Fig. |21 we show the lowest en- 
ergy structure along with some low lying excited state 
geometries of both clusters. The lowest energy geometry 
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FIG. 1: Specific-heats and average collision cross sections 
measured for size-selected Gair^ and Ga2o^ clusters as a 
function of temperature. The solid blue points show the 
specific-heats which are normalized to 3NkB (the classi- 
cal value), where fcs is the Boltzmann constant and A'^ = 
(3n — 6 -|- 3/2)/3 (n = number of atoms in the cluster, and 
3n — 6 and 3/2 are due to the vibrational and rotational con- 
tributions, respectively). 
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FIG. 2: The ground-state and some representative low-lying 
and excited-state geometries of Gair"*" and Ga2o'*^ clusters. 
Fig (1) corresponds to the ground-state geometry. The energy 
difi'erence AE is given in eV with respect to the ground-state. 
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of the Gaiv"*" cluster (see Fig. |21-a(l)) is similar to that 
of Gai7 reported in our earlier workJS It has a distorted 
decahedral structure, which suggests the possibility of 
further cluster growth to a 19~atom double decahedron. 
In contrast, the ground-state geometry of Ga2o^, shown 
in Fig. |21-b(l), is more symmetric. It can be described 
as a double decahedral structure of 19 atoms, with the 
bottom capped atom merging into the pentagonal plane 
to form a hexagonal ring. In addition, an atom from the 
top pentagon and the upper capped atom rearrange to 
accommodate the 20**^ atom, leading to a dome-shaped 
hexagonal ring. 

We now analyze the structural properties in detail to 
get an insight into the features that influence the melting 
characteristics. An analysis of the bond-length distribu- 
tion shows that there are 12 bonds, for each cluster, hav- 
ing distances less than 2.55 AiSl Interestingly, for Gai7+, 
these short bonds are spread all over the cluster, whereas 
for Ga2o^, they form the upper and the lower hexagonal 
rings. The distribution of coordination numbers^ indi- 
cate that for Ga2o^, almost all the atoms in the rings 
(about 16), have a coordination number of 4. The Gai7+ 
cluster, however, does not have such a uniform distribu- 
tion of coordination numbers. Thus, the ground-state ge- 
ometry of Ga.n'^ might be considered to be "disordered", 
while that of Ga2o^ exhibits a more-ordered structure. 

Striking differences are also observed in the low energy 
isomers and their distribution on the potential-energy 
surface. As mentioned above, we have obtained more 
than 50 distinct isomers spanning an energy range of 
about 1.0 eV above the ground-states for each cluster. 
In Fig. 121 we plot the energies of these isomers relative to 
the ground-state, arranged in an ascending order. The 
isomers for the Gai7"^ cluster appear to exhibit an almost 
continuous energy distribution. While a few of these iso- 
mers are severe distortions of the ground-state geome- 
try, the rest do not show any resemblance (see Fig. 
It appears that for the Gai7+ isomers in this low en- 
ergy regime, small rearrangements of the atoms, costing 
just a small amount of energy, lead to several close-lying 
isomers, so that the isomer distribution is almost con- 
tinuous. In contrast, the isomers of Ga2o"^ cluster are 
distributed in three groups, separated by an energy gap 
of about 0.2 eV (Fig. |3J). The first group of isomers have 
slightly different orientations of atoms in the hexagonal 
rings and are nearly degenerate with the ground-state. 
The second group consists of structures having only the 
lower hexagonal ring while the third group has no rings. 
This indicates that the hexagonal units of Ga2o^ clus- 
ter are stable and difhcult to break. The stability of 
the ring-pattern of Ga2o^ and the isomer distribution 
for both clusters should have a substantial effect on the 
melting characteristics. Indeed, as we shall see further 
below, these features play a crucial role in the finite- 
temperature characteristics. It should be mentioned that 
although these observations are based on rather limited 
search, we believe that the general features described here 
are essentially correct. 
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FIG. 3: The energies of the isomeric structures of Gai7^ and 
Ga2o'*' with respect to their ground-states. 



TABLE I: The number of basins with more than one atom 
at different values of the electron localization function for the 
ground-state structures of Gai7"'" and Ga2o^ clusters. The 
numbers in parenthesis represent the number of atoms in each 
basin. 



ELF value 


Gai7+ 


Ga2o''' 


0.85 





1 (2) 


0.77 


1 (2) 


2 (5,7) 


0.75 


3 (2,2,2) 


2 (5,7) 


0.73 


2 (3,4) 


1 (14) 



The most important difference between the two clus- 
ters is the nature of the bonding. We use the concept of 
an electron localization function (ELF),^^ to describe the 
nature of bonding. This function is normalized to a value 
between zero and unity; a value of 1 represents a perfect 
localization of the valence charge while the value for the 
uniform electron gas is 1/2. The locations of maxima of 
this function are called attractors, since other points in 
space can be connected to them by paths of maximum 
gradient. The set of all such points in space that are 
attracted by a maximum is defined to be the hasin of 
that attractor. Basin formations are usually observed as 
the value of the ELF is lowered from its maximum, at 
which there are as many basins as the number of atoms 
in the system. Typically, the existence of an isosurface 
or a basin in the bonding region between two atoms at a 
high ELF value, say > 0.70, signifies a localized bond in 
that region. 

We have analyzed the electron localization functions 
for Gai7"'" and Ga2o~'' clusters for values < 0.85. In Ta- 
bleQl we give the number of basins containing two or more 
atoms, for selected ELF values. The table clearly shows a 
fragmented growth pattern of the basins for Gai7"'", each 
containing very few atoms as compared to that of Ga2o~'". 
For instance, at an isovalue of 0.75, while Gai7"'' has three 



5 




correspond to merged basin structures. 



basins each having 2 atoms, Ga2o~'' has just two basins 
each containing 5 and 7 atoms that corresponds to the 
two hexagonal rings. The ELF contours for the isovalue 
of 0.75 are shown in Fig.^ The merged basins structures 
are shown by the black lines. It may be inferred that the 
bonds between atoms in the hexagonal rings of Ga2o"'' 
are strong and covalent in nature with similar strengths, 
while the fragmented basin growth pattern in Gary"*" in- 
dicates inhomogeneity of the bond strengths. 

The calculated, normalized, canonical specific-heats 
are shown in Fig. [S] plotted against temperature. The 
plot for Gary"*" exhibits a broad feature (apparently con- 
sisting of three components) which extends from 300 K 
to 1400 K. For Ga2o~'', the calculated specific-heat re- 
mains nearly flat up to about 600 K, it then increases 
sharply and peaks at about 800 K, in excellent agreement 
with the experimental results described above. Thus, in- 
teresting size-sensitive features seen in the experimental 
heat-capacities are reproduced in our simulations. This 
behavior can be understood from our earlier discussion 
of the bond-length distributions, coordination numbers, 
isomer-distributions, and the nature of bonding in these 
clusters. While the Ga.n'^ cluster shows no real evidence 
for ordered behavior, the Ga2o''' cluster has well-ordered 
ring-patterns. Thus, when Gair"*" is heated, the bonds 
soften gradually, and the cluster hops through all its iso- 
mers continuously. This is clearly demonstrated by the 
ionic motion as a function of temperature, which shows 
that this cluster evolves through all isomers smoothly 
from 300 K to 1400 K. On the other hand, the ionic mo- 
tion for Ga2o^ shows only minor rearrangements of the 
atoms until 600 K, and then the cluster visits all the iso- 
mers corresponding to the first group of isomers described 
above. At about 700 K, the upper hexagonal ring breaks, 
while at about 800 K, the lower ring breaks. Thus, melt- 
ing of Ga2o^ cluster is associated with the breaking of 
the well-ordered covalently bonded hexagonal units. 

We have also analyzed the melting characteristics via 
traditional parameters such as, the root-mean-squared 
bond-length-fluctuations ((5rms) and the mean-squared 
ionic displacements (MSD). In Fig. El we show the Srms 
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FIG. 5: Normalized canonical specific-heat for Gaiy"*" 
(continuous-line) and Ga2o^ (dashed-line). Co ~ (3n — 6-1- 
3/2)kB is the zero temperature classical limit of the rotational 
plus vibrational canonical specific-heat. 
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FIG. 6: Root-mean-square bond-length-fiuctuations (5rms) 
for Gai7^ (continuous-line) and Ga2o (dashed-line). 



for Gai7+ and Ga2o^ clusters. This plot correlates well 
with the specific-heat curve shown in Fig. [S] The S-cms 
for Gai7+ rises gradually from 300 K, while for Ga2o^, it 
rises sharply at about 700 K, and finally saturates to the 
same value for both clusters. It may be inferred from this 
observation that the behavior of both clusters at temper- 
atures say, T > 800 K, are similar and that both clusters 
can be considered to be in liquidlike states. This conclu- 
sion is further substantiated by the MSD plots (figures 
not shown), which saturate at sa 21A^, at about 1200 K 
for both clusters. 
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IV. SUMMARY AND CONCLUSION 

It is evident from the present study that the nature 
of the ground-state geometry and bonding strongly in- 
fluences the finite-temperature characteristics of Gai7+ 
and Ga2o'''- At high temperatures, T > 800 K, 
both Gai7+ and Ga2o~'' have similar root-mean-squared 
bond-length-fluctuations and the mean-squared ionic 
displacements so that both of them can be considered 
to be in liquidlike states. The experimental results show 
that while Gai7~^ apparently undergoes a solid-liquid 
transition without a significant peak in the speciflc-heat, 
Ga2o^ melts with a relatively sharp peak. The simula- 
tions show that if the cluster is "ordered" (i.e. a large 
fraction of the constituent atoms show similar bonding, 
coordination numbers, and bond energies) then it is likely 
to show a sharp melting transition with a significant peak 
in the speciflc-heat. On the other hand, if the cluster is 
"disordered" (i.e. the constituent atoms occur in a wide 
distribution of bonding environments) it will probably 
undergo a solid-liquid transition without a significant 
peak in the speciflc-heat. In the latter case, the num- 
ber of isomers or conformations sampled by the cluster 
increases steadily as the temperature is raised, instead of 
the abrupt increase that occurs when a cluster undergoes 
a sharp melting transition. 

These observations have interesting consequences for 
the flnite-temperature behavior of small clusters as a 
function of cluster size. It is likely that as clusters grow in 
size their structures evolve from one well-ordered struc- 



ture to another, passing on the way through some clus- 
ter sizes that have "disordered" structures. For instance, 
the 13-atom gallium cluster is predicted to have a highly 
symmetric decahedron structure with a bonding pattern 
that is similar to that found here for Ga2o^."'^^ So in the 
present case, cluster growth from Gai3 (a decahedron) 
to (Ga2o'^, a distorted double-decahedron) proceeds via 
a disordered Gaiy^ structure. Such behavior is also ob- 
served for sodium clusters in 40 to 55 atom size range; 
the ground-state geometries of Na4o and Nass are ei- 
ther icosahedron or close to icosahedron while that of 
Na5o has no particular symmetryii2ii& In such cases, we 
expect that the speciflc-heats should change from show- 
ing a well-defined peak to a rather broad one, and back 
again to well-defined. We believe this behavior to be 
generic as it has not only been observed in case of gal- 
lium clusters*'^ but also in case of aluminum clusters^fi, 
experimentally, and for sodium clusters in the simula- 
tions mentioned above. 
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